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The  spectral  density  S^{f)  =  A^liJ  jl  Hz)“  of  magnetic  flux  noise  in  ten  dc  superconducting  quantum 
interference  devices  (SQUIDs)  with  systematically  varied  geometries  shows  that  a  increases  as  the 
temperature  is  lowered;  in  so  doing,  each  spectrum  pivots  about  a  nearly  constant  frequency.  The  mean- 
square  flux  noise,  inferred  by  integrating  the  power  spectra,  grows  rapidly  with  temperature  and  at  a  given 
temperature  is  approximately  independent  of  the  outer  dimension  of  a  given  SQUID.  These  results  are 
incompatible  with  a  model  based  on  the  random  reversal  of  independent,  surface  spins. 

DOI:  10.1103/PhysRevLett.ll0.147002  PACS  numbers:  85.25.Dq,  05.40.-a,  85.25.Am 


Flicker  (1//)  noise  at  low  frequencies  /  is  observed  in 
all  solid  state  devices.  In  particular,  1//  noise  in  super¬ 
conducting  quantum  bits  (qubits)  [1]  is  the  leading  cause  of 
decoherence  (loss  of  phase  coherence)  as  the  state  of  the 
qubit  evolves.  The  resulting  reduction  in  the  pure  dephas¬ 
ing  time  T0  limits  the  fidelity  of  qubit  readout  and  inhibits 
scaling  up  to  multiple  qubit  circuits.  In  some  cases,  the 
origin  of  the  noise  is  reasonably  well  understood.  For 
example,  in  charge-sensitive  devices  such  as  charge  qubits 
[2]  1//  noise  arises  from  the  hopping  of  electrons  between 
nearby  traps  that  induce  fluctuating  charges  on  the  qubit.  In 
the  case  of  Josephson  tunnel  junctions — which  are  constit¬ 
uents  of  flux  [3]  and  phase  [4]  qubits  and  superconducting 
quantum  interference  devices  (SQUIDs)  [5] — 1//  noise  in 
the  critical  current  (maximum  supercurrent)  arises  from  the 
trapping  of  electrons  in  the  barrier  and  their  subsequent 
release  [6].  The  mechanism  for  magnetic  flux  1//  noise, 
however,  with  a  power  spectrum  =  AV(//1  Hz)“, 
has  remained  poorly  understood.  Here,  O  is  flux  and 
a  £  1 .  Flux  noise  is  a  major  source  of  intrinsic  dephasing 
[7]  in  superconducting  flux  [8-12]  and  phase  qubits 
[7,13,14]  and  in  the  quantronium  [15];  it  also  generates 
low-frequency  noise  in  SQUIDs  [16-18]. 

It  has  been  proposed  that  flux  noise  originates  from  the 
random  reversal  of  spins  at  the  interfaces  between  the 
superconductor  and  its  substrate  and  surface  oxide  layer 
[19-21].  Models  that  assume  independent  spins  with  a 
magnetic  moment  /Zg  (the  Bohr  magneton)  predict  an 
areal  spin  density  of  5  X  10*^  m“^.  The  same  value  is 
found  from  measurements  of  paramagnetism  in  SQUIDs 
[22]  and  Au  rings  [23].  Noise  correlation  measurements 
on  flux  qubits  are  consistent  with  the  surface  spin 
model  [11,12]. 

A  crucial  question  is  whether  the  independent-spin 
model  describes  the  flux  noise  dynamics  correctly  or  if 


spin-spin  interactions  are  significant.  Sendelbach  et  al.  [24] 
inferred  the  existence  of  clusters  from  measurements  of 
inductance  fluctuations  in  SQUIDs,  suggesting  a  non- 
negligible  spin-spin  interaction.  The  independent-spin 
model  was  used  in  an  analytical  calculation  of  the  mean- 
square  value  of  the  flux  noise  (O^)  as  a  function  of  loop 
geometry  [13,20].  Measurements  of  decoherence  in  flux 
qubits  [10]  were  in  agreement  with  the  predicted  scaling 
while  measurements  of  decoherence  in  phase  qubits 
instead  showed  a  scaling  with  loop  inductance  [14]. 

In  this  Letter  we  report  flux  noise  measurements  as  a 
function  of  temperature  for  ten  dc  SQUIDs  with  system¬ 
atically  varied  geometries.  We  find  that  both  Stjjjl  Hz)  and 
the  slope  a  vary  systematically  with  temperature. 
Remarkably,  for  a  given  SQUID  S^{f)  pivots  about  a 
nearly  fixed  frequency  as  the  temperature  is  changed. 
Values  of  (O^)  inferred  from  our  measurements  of  S$(/) 
deviate  markedly  from  the  predicted  scaling  with  loop 
geometry.  Furthermore,  although  there  is  no  evident  tem¬ 
perature  dependence  in  the  prediction  of  over  the 
range  of  temperatures  investigated,  pivoting  causes  the 
inferred  values  of  (O^)  to  vary  over  several  orders  of 
magnitude  as  we  change  the  temperature.  We  are  unable 
to  reconcile  our  data  with  the  independent-spin  model. 

We  measured  S,p{f)  in  a  total  of  20  dc  SQUIDs,  con¬ 
nected  in  series  in  sets  of  5  on  four  chips,  fabricated  on 
oxidized  Si  wafers  using  a  trilayer  Nb-AlOx-Nb  junction 
technology.  Two  chips  were  made  at  MIT-LL  and  two  at 
NIST.  Since  the  two  chips  from  each  institution  showed 
very  similar  results,  we  present  data  on  only  one  from  each, 
labeled  I  (LL)  and  II  (NIST).  The  SQUIDs  were  patterned 
in  a  square  geometry  with  outer  widths  2R  and  loop  line 
widths  W  listed  in  Table  I,  along  with  their  estimated  loop 
inductances.  Each  junction  was  resistively  shunted  to 
eliminate  hysteresis  on  its  current-voltage  characteristic. 
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TABLE  I.  Dimensions  and  inductances  of  SQUIDs  I  and  11. 


R  (/um) 

IT  (/rm) 

R/W 

L  (pH) 

I.l 

12 

0.5 

24 

80 

1.2 

6 

0.5 

12 

33 

1.3 

3 

0.5 

6 

12 

1.4 

1.5 

0.5 

3 

4 

1.5 

1.5 

0.5 

3 

4 

II.  1 

265 

240 

1.1 

120 

II.2 

145 

120 

1.2 

98 

II.3 

85 

60 

1.3 

92 

II.4 

55 

30 

1.8 

96 

II.5 

40 

15 

2.7 

106 

As  shown  schematically  in  Eig.  1,  the  hve  SQUIDs  on  a 
given  chip  were  connected  in  series  with  an  off-chip  com¬ 
pensating  resistor  =  0.46  H,  the  input  coil  of  a  readout 
SQUID  operated  in  a  flux-locked  loop  (ELL),  and  a  choke 
to  suppress  high-frequency  currents  [25].  The  entire  circuit 
was  enclosed  in  two  superconducting  shields  and  a  cryo- 
perm  shield  to  reduce  ambient  magnetic  field  fluctuations 
and  the  Earth’s  static  magnetic  field.  The  magnetic  flux  in 
the  SQUIDs  was  established  by  means  of  a  current  in  a 
common  coil.  To  measure  the  flux  noise  in  a  given  SQUID, 
we  biased  it  at  a  voltage  of  2.5  or  5  /U,V  with  a  current  7^ 
and  canceled  the  quasistatic  current  induced  in  the  input 
coil  and  remaining  SQUIDs  with  a  compensating  current 
in  R^.  Fluctuations  50  in  O  of  the  measured  SQUID 
induced  a  current  in  the  input  coil  of  the 

readout  SQUID.  For  all  data  represented  here,  we  biased 
the  SQUID  where  dl/ dO — determined  by  measuring  the 
response  to  a  small  oscillating  flux — was  a  maximum.  The 
temperature  T,  measured  using  calibrated  Ge  (mixing 
chamber)  and  RuOx  (sample  box)  resistance  thermome¬ 
ters,  was  stabilized  with  feedback  from  the  Ge  thermome¬ 
ter  to  better  than  1  part  in  10^  during  data  acquisition  [26]. 
Nyquist  noise  from  R^  yielded  a  temperature  within  ±5% 
of  that  of  the  thermometers. 

We  acquired  a  time  series  of  the  voltage  fluctuations  for 
1  hour,  computed  the  spectral  density  and  converted  it  to 


FIG.  1  (color  online).  Configuration  of  measurement  system  to 
measure  flux  noise  in  three  SQUIDs.  The  bias  current  If,  enables 
the  measurement  of  the  middle  SQUID.  The  static  voltage  across 
the  SQUID  is  canceled  by  the  current  R  applied  to  R^. 


>5<i>(/)-  We  performed  a  least  squares  fit  to  S^(f)  = 
Hz)"  -f  C^,  representing  the  flux  1//"  noise 
and  the  white  noise  from  the  resistive  shunts,  to  obtain  A 
and  a.  To  confirm  that  the  inferred  value  of  A^  was 
independent  of  dljd^,  we  measured  the  noise  at  other 
values  of  O  and  dl/d^.  In  addition,  we  made  measure¬ 
ments  at  dlj d^  =  0,  enabling  us  to  determine  the  critical 
current  1//  noise.  We  verified  that  the  power  spectrum  of 
the  white  noise  from  the  resistive  shunts  dominated  those 
from  R^  and  the  readout  electronics,  and  that,  with  a  flux 
bias  of  nOo,  it  remained  white  at  frequencies  down  to  the 
1//  knee  of  the  readout  electronics  with  a  magnitude  in 
reasonable  agreement  with  theoretical  predictions  [27]. 
Here,  n  is  an  integer  and  Oq  =  /i/2e  is  the  superconduct¬ 
ing  flux  quantum. 

As  representative  data,  in  Fig.  2(a)  we  show  raw  power 
spectra  for  SQUID  II.3.  The  spectra  are  1  / /-like,  with  a 
slope  that  flattens  at  higher  frequencies  as  the  white 
noise  from  the  shunt  resistors  becomes  significant.  At 
low  frequencies  (/  £  10“'  Hz)  and  high  temperatures 
(T  S  1.2  K),  fluctuations  in  the  critical  current  are  signifi¬ 
cant,  thereby  increasing  the  slope  of  the  spectra  [26]. 
Consequently,  in  our  fits  to  the  spectra,  we  disregard  noise 
from  this  region. 

The  fitted  values  of  P?  and  a  are  plotted  in  Fig.  3  for  all 
five  SQUIDs  on  both  chips.  We  see  immediately  that  the 
systematic  trends  with  T  for  each  chip  differ  between  the 
two  chips.  For  chip  I,  A^  increases  as  T  is  lowered  to  0. 1  K, 


c^O 

le 


CO 


FIG.  2  (color  online).  Power  spectra  S^(f).  (a)  As-acquired 
flux  noise  spectral  densities  for  SQUID  II.3  at  1 1  temperatures. 
Inset  shows  data  at  0.1  K  with  corresponding  fit  to 
4^/(//l  Hz)"  -I-  (dashed  line).  Fits  of  S^)/)  for 
(b)  SQUID  I.l  at  10  temperatures  and  (c)  SQUID  II. 3  at  11 
temperatures. 
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FIG.  3  (color  online).  Temperature  dependence  of  1//  flux 
noise  for  10  SQUIDs.  (a),  (b)  Fit  coefficients  for  SQUIDs  I 
and  II  vs  T.  (c),  (d)  Fit  coefficients  a  for  SQUIDs  I  and  II  vs  T. 
Confidence  in  fits  is  ±10%  in  (a),  (b);  ±0.03  (T  <  1.1  K)  and 
±0.05  (7’>  1.1  K)  in  (c);  and  ±0.03  (T  <  1.1  K)  and  ±0.07 
(r>  1.1  K)  in  (d). 

with  an  upturn  for  T  <0.15  K,  and  increases  monotoni- 
cally  with  R.  In  contrast,  for  chip  II hecomes  indepen¬ 
dent  of  T  for  r  £  0.6  K,  a  result  reminiscent  of  the 
findings  of  Wellstood  et  al.  [18].  With  regard  to  the  slope, 
for  both  chips  a  increases  dramatically  as  T  is  lowered 
from  4.0  to  0. 1  K.  For  chip  I,  the  value  of  a  increases  with 
R.  Similarly  low  values  of  a,  0.4-0.5  at  4  K,  have  been 
observed  by  other  authors  [17,28].  Although  we  have  no 
explanation  for  the  disparate  dependencies  of  and  a  on 
T,  we  remark  that  for  chip  I  VP  is  constant  at  0.5  /xm  while 
for  chip  II  W  varies  from  15  to  240  /xm.  A  conceivable 
explanation  is  that  interactions  between  spins  give  rise  to 
spatial  correlations. 

The  progressive  increase  in  a  as  T  is  lowered  is  illus¬ 
trated  vividly  in  Figs.  2(b)  and  2(c),  where  we  plot 
AV(//1  Hz)“  for  SQUID  1.1  [Fig.  2(b)]  and  11.3 
[Fig.  2(c)]  for  all  T  at  which  we  obtained  data.  Ignoring 
the  two  highest  spectra  in  Fig.  2(b) — corresponding  to 
7’<  0. 15  K  with  upturns  in  the  data  in  Fig.  3(a) — quite 
unexpectedly  and  remarkably  we  observe  that  each  set  of 
spectra  pivots  around  an  approximately  constant  frequency 
fp\  fp  =  16±  20  Hz  in  Fig.  2(a)  and  fp  =  6.2  ±  0.5  Hz 
in  Fig.  2(b).  The  Supplemental  Material  [29]  shows  similar 
plots  for  all  ten  SQUIDs.  The  plots  in  Fig.  2  imply  that, 
although  one  traditionally  defines  the  noise  magnitude  as 
A^  =  5'<i,(l  Hz),  in  fact  A^  is  an  incomplete  characteriza¬ 
tion  in  the  absence  of  a  knowledge  of  a. 

We  next  estimate  the  mean-square  flux  noise  (O^).  For  a 
circular  loop  of  outer  radius  R  and  line  width  W,  in  the 
limit  R  »  W,  Bialczak  et  al.  [13]  showed  that 

2al  ,  R  V\n{2bW/\^)  1 


for  uncorrelated  spins  with  uniform  surface  density  a  and 
magnetic  moment  /x^.  Here,  /xq  is  the  vacuum  permeabil¬ 
ity,  b  is  the  film  thickness,  and  A  is  the  penetration  depth. 
For  the  square  geometry  of  our  devices,  there  is  a  multi¬ 
plicative  correction  of  order  unity,  which  we  neglect. 
Because  we  cannot  directly  measure  (O^),  we  use 

(02)=  (2) 
Jfl 

to  relate  (O^)  to  our  measurements  of  Here,  we 

extrapolate  S^{f),  which  we  measure  typically  over  sev¬ 
eral  decades,  to  /i  and  /2.  We  set  fi  =  10“"^  Hz,  approxi¬ 
mately  the  lowest  frequency  to  which  flux  noise  has  been 
measured.  The  value  of  /2  is  more  difficult  to  estimate 
as  it  is  poorly  understood.  Since  Bylander  et  al.  [30]  and 
Slichter  et  al.  [31]  observed  flux  noise  up  to  measurement- 
limited  frequencies  of  20  MHz  and  1  GHz,  respectively,  we 
use  /2  =  10^  Hz  for  our  analysis.  Separate  measurements 
on  Al  SQUIDs,  fabricated  in  exactly  the  same  way  in  the 
same  equipment  as  the  qubits  in  Ref.  [31],  yielded  average 
low-temperature  values  remarkably  similar  to  those 
observed  in  SQUIDs  I  and  II:  A^  =  2. l(/xOo)^/Hz  and 
a  =  0.72.  These  results  suggest  that  details  of  the  process¬ 
ing  do  not  greatly  alter  the  characteristics  of  flux  noise.  In 
particular,  we  expect  /2  to  be  comparable  in  SQUIDs  and 
qubits  alike. 

For  all  ten  SQUIDs,  Fig.  4  shows  inferred  using 

Eq.  (2)  and  the  data  plotted  in  Fig.  3,  vs  T.  For  comparison, 
the  dashed  lines  indicate  the  value  of  for  each 

SQUID  calculated  using  Eq.  (1),  the  values  of  R  and  W 
from  Table  I,  a  =  5  X  10'^  m“^,  A  =  39  nm,  and  b  = 
150  nm  (I),  200  nm  (II).  While  is  independent  of 

r  for  r  T^,  the  values  of  increase  strongly  with 

increasing  T.  At  the  lowest  T,  values  of  exceed  those 
of  by  1  to  2  orders  of  magnitude  and  at  the  highest  T 

by  3  to  5  orders  of  magnitude.  Clearly,  the  rapid  calculated 
increase  of  with  T  is  inevitable  given  the  behavior 

in  Fig.  2. 

As  a  further  test  of  the  theory,  we  investigated  the 
scaling  of  S^{f)  with  SQUID  dimensions  to  compare 


FIG.  4  (color  online).  Mean-square  flux  noise,  inferred  from 
measured  power  spectra  using  Eq.  (2),  vs  T.  (a)  Chip  I  and 
(b)  Chip  II.  Horizontal  dashed  lines  represent  the  predictions  of 
Eqs.  (1)  and  (2)  for  the  values  of  R  and  W  listed  in  Table  I. 
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with  the  predictions  of  Eq.  (1).  For  brevity,  we  present 
results  only  for  SQUIDs  I,  for  which  W  is  constant,  so  that 
(cp^aic)  Figurs  5(a)  shows  5<i,(54  Hz)  vs  R  for  the  four 
distinct  geometries  of  SQUIDs  1.  Here,  we  chose  the 
average  value  of  54  Hz,  to  minimize  the  effect  of  the 
T  dependence  of  the  spectra  on  the  geometric  scaling.  To  a 
reasonable  approximation,  apart  from  the  data  for  SQUIDs 
1.1  and  1.2  at  the  two  lowest  temperatures,  5'$(54  Hz) 
scales  linearly  with  R.  Figure  5(b)  shows  vs  R  for 

ten  temperatures,  which  we  see  is  approximately  indepen¬ 
dent  of  R,  increases  with  T,  and  has  a  magnitude  substan¬ 
tially  higher  than  the  prediction. 

We  now  discuss  potential  sources  of  error  in  that 

would  give  rise  to  a  strong  temperature  dependence.  A 
T-dependent  spin  density  is  precluded  by  the  Curie-Law 
behavior  of  the  spin  paramagnetism  [22,23].  We  also  note 
that,  for  the  observed  values  of  a,  {‘Ffnf)  changes  by  less 
than  1%  as  /i  is  varied  from  zero  to  10“^  Hz.  With  regard 
to  /2,  one  could  argue  that  the  value  of  /2  =  1  GHz, 
observed  at  50  mK  [31],  becomes  lower  as  T  increases. 
However,  to  reduce  the  highest-temperature  value  of  (O?  f) 
in  Fig.  5  to  its  lowest  temperature  value  by  lowering  /2 
would  require  /2  to  be  roughly  700  kHz  for  SQUIDs  I  and 
just  8  kHz  for  SQUIDs  11.  Such  values  of  /2  are  improbably 
low;  it  is  difficult  to  imagine  a  noise  process  that  extends  to 
1  GHz  at  0. 1  K,  yet  to  less  than  1  MHz  at  4  K.  Regarding  a, 
high-frequency  measurements  in  qubits  have  sampled  the 
flux  noise  spectrum  from  0.2  to  20  MHz  [30]  as  well  as 
from  10“2  to  1  Hz,  1  to  20  MHz,  and  0.7  to  1  GHz  [31]. 
These  measurements  each  align  with  a  spectrum  of  con¬ 
stant  a  {a  =  0.9  [30]  and  a  =  0.57  [31])  and  value  of 
54,(1  Hz)  typical  of  those  observed  in  SQUIDs.  These 
measurements  suggest  that,  while  a  can  vary  between 
devices,  within  a  single  device  it  is  frequency  independent 
over  a  very  wide  bandwidth.  Finally,  errors  in  our  fits  of  a 


R  (//m) 

FIG.  5  (color  online).  Scaling  of  noise  power  with  R  for 
SQUIDs  I  for  ten  temperatures  from  0.1  K  to  4.0  K. 
(a)  54,(54  Hz)  vs  R.  Dashed  line  indicates  scaling  with  R,  with 
arbitrary  magnitude,  (b)  vs  R.  Dashed  line  represents 

prediction  of  Eq.  (1). 


will  inevitably  lead  to  errors  in  For  instance,  at  low 

(high)  temperatures  decreasing  a  by  0.03  (0.07)  increases 
)  by  a  factor  of  roughly  1 .6  (4).  These  errors,  while  not 
negligible,  cannot  account  for  the  strong  temperature  de¬ 
pendence  and  discrepancies  of  several  orders  of 
magnitude. 

Given  the  difficulty  of  reconciling  our  data  with  the 
predictions  of  a  model  based  on  N  single,  uncorrelated 
spins — for  which  (O^)  (x  A/x| — we  consider  the  possibil¬ 
ity  that  the  spins  form  clusters  [24].  We  assume  the  clusters 
to  be  uncorrelated  with  each  other  and  to  contain  an 
average  of  Z  spins  producing  a  magnetic  moment  /x^, 
where  Z  may  depend  on  T.  There  are  three  scenarios: 
(i)  ferromagnetic  clusters  with  /x^  =  Z/Xg  and  (O^)  0^ 
(A/Z)(Z/(Xg)^  =  AZ/x|,  (ii)  random  (glassy)  clusters 
with  /x^  =  Z'/^B  and  (O^)  a  (A1/Z)(Z‘/2/xb)^  = 
and  (iii)  antiferromagnetic  clusters  with  /x^,  ~  /x^  and 
(02^  ot;  (^/2)/x|.  As  Z  grows,  the  rate  at  which  a  cluster 
reverses  its  magnetic  moment  is  assumed  to  decrease 
rapidly  [32].  At  a  given  T,  such  a  distribution  of  reversal 
rates  offers  a  natural  explanation  for  the  13-decade  range  in 
lifetimes  required  for  S^if)  to  range  from  10“^^  to  10®  Hz 
[33].  Furthermore,  a  plausible  scenario  for  the  spectral 
pivoting  as  T  is  lowered  is  a  progressive  growth  in  Z, 
reducing  the  number  of  small  clusters  generating  noise 
above  fp  and  increasing  the  number  of  larger  clusters 
generating  noise  below  fp — thereby  increasing  a.  This 
simplistic  picture  has  implications  for  the  temperature 
dependence  of  each  of  the  three  cases.  As  T  is  lowered, 
(O^)  increases  for  ferromagnetic  clusters,  remains  constant 
for  random  clusters,  and  decreases  for  antiferromagnetic 
clusters. 

Finally,  the  noise  behavior  must  be  compatible  with  the 
observed  spin  paramagnetism  [22,23],  which  follows  a 
Curie  \/T  scaling.  We  note  that  the  contribution  of  a 
cluster  with  magnetic  moment  /x^.  to  the  linear  magnetic 
susceptibility  scales  as  On  the  other  hand,  once 

fji^B  >  kgT  the  susceptibility  ceases  to  be  linear  in  mag¬ 
netic  field,  the  magnetization  of  the  cluster  saturates  and 
Curie  scaling  is  violated.  Consequently,  since  Curie  behav¬ 
ior  was  observed  at  T  =  50  mK  and  5=10  mT,  we  find 
/x^  <  7/Xb.  Thus,  Z  <  7  for  the  ferromagnetic  case,  Z  <  7^ 
for  the  random  case,  and  Z  is  unrestricted  for  the  antifer¬ 
romagnetic  case.  We  also  point  out  that  the  areal  single 
spin  density  of  5  X  10'^  m“^  derived  from  both  noise  and 
paramagnetism  experiments  is  unchanged  for  glassy  clus¬ 
ters.  Needless  to  say,  a  detailed  understanding  of  the  inter¬ 
actions  between  both  the  spins  that  form  such  clusters  and 
between  the  clusters  themselves  will  be  required  to  develop 
a  credible  explanation  of  flux  noise. 
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